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Abstract 

This note deals with estimating the volume of the set of separable mixed quantum states when 
the dimension of the state space grows to infinity. This has been studied recently for qubits; here we 
consider larger partides and conclude that, in all cases, the proportion of the states that are separable 
is super-exponentially small in the dimension of the set. We also show that the partial transpose 
criterion becomes imprecise when the dimension increases, and that the lower bound Q~ N / 2 on the 
(Hilbert-Schmidt) inradius of the set of separable states on N qubits obtained recently by Gurvits 
and Barnum is essentially optimal. We employ Standard tools of classical convexity, high-dimensional 
probability and geometry of Banach spaces. One relatively non-standard point is a formal introduc- 
tion of the concept of projective tensor products of convex bodies, and an initial study of this concept. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.65.Db, 02.40.Ft, 02.50.Cw 

1 Introduction and summary of results 

An important problem in quantum information theory is to estimate quantitatively parameters related 
to entanglement. This phcnomenon is thought to be at the heart of quantum information processing 
whilc, on the other hand, its experimental creation and handling are still challcnging. The question of 
determining the volume of the set of separable (i.e., uncntanglcd) states has been asked for example in 
^12], where we refer the reader for background and motivation. In this paper we obtain estimates which 
are meaningful when the dimension of the state space is "not too small," and which are asymptotically 
tight as that dimension tends to infinity. The effcctive radius in the sense of volume of (that is, the radius 
of the Euclidean ball of the same volume as; also referred to as volume radius) the set of separable states 
for qubits has been determined (prcciscly on the scale of powers of the dimension) in Pj. Here we shall 
deal with more general partides, while using the same general approach via the methods of asymptotic 
convex geometry. The conclusion is that when the dimension of the state space is large, then all but 
extremely few (as measured by the Standard volume) states are entangled; see Theorems 1 and 2 below 
for precise statements. As further application of our techniques we show the optimality (up to factors of 
lower order) of the lower bound on the radius of a Hilbert-Schmidt ball containcd in the set of separable 
states on N qubits that was obtained recently by Gurvits and Barnum 0]. Still another consequence is 
that, for large dimensions, the set of states with positive partial transpose jS] |2] is, in the sense of the 
volume radius, much closer to the set of all states than to the set of separable states; in other words, the 
Peres partial transpose criterion becomes less and less precise as the dimension increases. 

We now recali the mathematical framework and introduce some notation. The state space is a complex 
Hilbert space H = C Dl <g> ■ • • ® C° N . We writc d = D 1 . . . D N for the dimension of H. This Hilbert space 
allows to describe quantum interactions between N partides; Dj = 2 corresponds to qubits and Dj = 3 
to qutrits. We write B(H) for the space (or C*-algebra) of linear maps from Tí into itself, and £? sa (7í) for 
the (real linear) subspace of self-adjoint operators. The space B(1i.) is cndowed with the Hilbert-Schmidt, 
or Frobenius, sealar product defined by (A, B)hs : = trA^B ( = Ylij=i AijBij if A and B are represented 
as matrices). A state on B(TÍ) (which we will abbreviate to "a state on H") can be represented as a 
positive (semi-definite) trace one element of B sa {TL). A state is said to be separable if it can be written 
as a convex combination of tensor products of N states, otherwise it is called entangled. That is, the set 



of states (also called in this context density operators or density matrices) is 

V = V{H) :={p£B sa {H),p^O,trp=l}, 

and the set of separable states is 

S = S{H) := conv{pi ® • • • ® p N) pj G P(C^)}. 

The notation S(7í) is in principio ambiguous: separability of a state on B(7í) is not an intrinsic property 
of the Hilbert space Ti or of the àlgebra B{7i)\ it depends on the particular decomposition of Tí as a 
tensor product of (smaller) Hilbert spaces. However, this will not be an issue here since our study focuses 
on fixed decompositions. 

T> and S are convex subsets of B sa (TL) of (real) dimension n := d 2 — 1. We write T\ for the affine 
subspace generated by T> (or S); it is the hyperplane of trace one matrices. The space B sa (Ji) inherits a 
(real) Euclidcan structure from the scalar product (-, -)hs- The corrcsponding unit ball will be denotcd 
by Bhs (or Bjj S ). We denote by vol the corrcsponding d 2 -dimcnsional Lebcsguc volume on B sa (TÍ); we 
will also write vol for the n-dimensional volume induced on T\ . 

We shall concentrate on the special case when all the subsystems are identical, i.e., Di = ■ ■ ■ = = 
D. Our techniques also apply to general state spaces, but formulae which fit all the situations would be 
more cumbersome than for this "homogeneous" case. If H = (C D )® N , there are two "regular" ways to 
make its dimension tend to infinity: either TV ^ 2 is fixed and D tends to infinity, or D ^ 2 is fixed and 
./V goes to infinity (as in [^j for D = 2). Surprisingly, it turns out that in the latter case the natural 
Euclidean structure is not the one given by the Hilbert-Schmidt scalar product, but rather by a mixture 
of that product with a Killing-type form. Our goal is to prové the following theorems (recali that we 
work with H = (C D )® N and that d = dímH = D N , TV = log D d and n = dimV = d 2 - 1 = D 2N - 1). 

Theorem 1 (Small number of large subsystems) There exist universal constants C, c > such 
that, for all D, N > 2, 

c N „ / vd1(5) \ 1/w C(NlnN) 1 / 2 

d l/2-l/2N < \^oï(V)J ^ dl/2-l/2JV 

We point out that considering the nth root of the ratio of volumes is natural from the geomètric point 
of view: if the two sets had been Euclidean balls, we would have obtained the ratio of their radii, and so 
we are comparing the volume radii of S and T> (see section ^.ll for additional geomètric background). 

Theorem 2 (Large number of small subsystems) There exist universal constants C',c' > such 
that, for all D,N ^ 2, 

d „ / vol(«S) \ l,n C'(DN\nN) 1 / 2 

d l/2+ aD ^ ^vol(O) ) ^ rf!/2+a D 

where 

a D ■= \ l0g D (l + -) - ^5 10 &d(D + 1) 

For illustration, we list approximate vàlues of «d for small D's: a-i « 0.094, as 0.061, « 0.044. Note 
that ao is a positive decreasing function of D, asymptotically equivalent to l/(2DlnD). The quantity 
d 1/2+OD may be alternatively expressed as ((D + l) 1 " 1 / d2 ) n / 2 . 

Theorem 1 asserts that, as D — > co, the ratio of the volume radii of S and T> is of order (f - ^ 1 / 2-1 / 2 ^ = 
jj-(n-i)/2 £ Q a mu itiplicative constant dcpcnding only on N). Theorem 2 is slightly less definitc: 
whilc, for fixed D, the ratio in question is determined prcciscly on the sealc of the powcrs of d, there are 
some (possibly parasitic) logarithmic factors. [We recali again that TV = log^, d — Ind/lnD.] We point 
out that, apart from the value of the numcrical constant, the upper estimate in Theorem 2 is always 
sharper than that in Theorem 1. However, the gain au in the exponent is negligible if N is fixed and 
D -> oo. 

It is possible to give reasonable estimates on the constants appearing in both theorems. One can, for 
example, take c = l/y/6, c' = 0.3 and C = C = 4.4. It is also possible to establish a tighter asymptotic 
behavior of these constants (relevant if one is only interested in large vàlues of d and particularly TV): 
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one can have d tending to e^^/yfh: m 0.845 as D or N tend to infinity, and C and C tending to 
V^e 1 / 4 ~ 1.816 as N tends to innnity. Finally, sharper estimates for specific dimensions can be obtained 
by using known results on, among others, efhcient spherical codes; see Appendix G of [3] for an example 
of such calculation. 

The case D = 2 of Theorem 2 was treated in detail in [2j . In this note we present additional ingredients 
required to deal with the general case. We point out that the (relatively) easy part is the estimation of 
vol(P), which of course does not depend on the tensor structure oiTÍ. It has been shown in [ü] that 

vol(P)^Vd(2 7 r)^- 1 )/ 2 r(1 ^ (d) (1) 

A direct calculation then yields that the volume radius of V behaves as e~ 1 / 4 £f _1 / 2 (l + 0(l/d)) as d — ► oo. 
It is also possible (and easy) to estimate the volume radius of T> using techniques with the same flavor as 
the ones presented here. For example, it was shown in [3] that for any value of d that radius is contained 
between gï _1 / 2 /2 and 2d~ x / 2 (see section 12721 below) . Finally, the fact that the volume radius in question 
is of order g? -1 / 2 as d — > oo follows from an early paper [7]- 

The ratio of Euclidean volumes is not the only geomètric parameter which is of interest here. An 
arguably more relevant parameter would be its analogue for the so-called Bures volume, which may be 
a more appropriate measure of size in the present context (see Another quantity with physical 

significance is the largest number e = e{TL) below which the mixturo (1 — e)Id/d + ep is separable for any 
state p. Of course it is possible to derive from volume estimates an upper bound on this e(Tl), but this 
bound is weaker that what is already known (|§1E3E])- Still, the knowledge of the volume radii yields 
additional information also in this context: for e ^> (vol( l S)/vol(2?)) 1 ^ dlmI ', the mixturo (1 — e)Id/d + ep 
is entangled for "most of p e V, and not just for very special states as those constructed in P ITTil 14]. 
Another quantity that has been studied in |1U1 21 in] IS the inradius of S, which is the radius of the 
inscribed Euclidean ball in S (here we assume that the balls are centered at Id/d, which is the natural 
center of both S and T>). It seems that, prior to results on qubits in the only upper bound which 
was known on the inradius of S was the inradius of D, or 1/ y/d{d — 1). [A remarkable fact discovered 
in jl·l| is that the two inradii coincide when N = 2.] Theorem 2 improves on this for large N giving, 
up to logarithmic factors, the upper bound l/d 1+a °. Moreover, the argument shows that the same is 
true for the symmetrized body E := conv(S U —S). [Cf. section l2~2l and note that the distance from 
the hypcrplane containing S to the origin is 1/yd > l/y/d(d — 1), so the inradius of E is at least as 
large as that of S.] The techniques developed in jTJj and in the present paper allow to tighten this bound 
substantially. We present the following sample result. 

Theorem 3 (Tight upper bounds on the inradii of S and E) There exists an absolute constant 
Co > such that, for any N, the set Y*((C 2 )® N ), and hence also S((C 2 )® N ), does not contain any 
Hilbert-Schmidt ball of radius C VN\tíN 6~ n / 2 . 

For comparison, lct us cite the lower bound of 6~ w / 2 on the inradius of the set E noted in Appendix 
H, and a similar lower bound on the a priori smaller inradius of S shown subscqucntly in 0]. This means 
that, at least for N qubits, the inradii of S and of E have been determined up to a multiplicative factor 
which is logarithmic in the dimension of the state space. The argument we present can be generalized 
to other vàlues of D, leading to an upper bound for the inradii of E((C D ) 8,Ar ) and <S((C· D )® AÍ ) which is 
(D(D + 1)) _JV / 2 (again, up to a logarithmic factor). For the inradius of the set E this bound is essentially 
sharp. However, for D > 2, the value (D(D + \))~ N / 2 does not coincide with the lower bound for the 
inradius of S((C 2 )® N ) that was found in g| and which is of the form 0((D(2D - 1)) _AÍ/2 ). As was 
noticed in 0], a better upper bound for D > 2 is derived from our upper bound for D = 2 and is (up 
to a logarithmic factor) of the form (1.5 x D 2 )~~ N / 2 . However, this still leaves a gap between the upper 
and lower estimates which is exponential in (as N — > oo). A similar discussion of the bilateral case 
(N = 2, D — ► oo) can be found in the Remark at the end of section 3. 

Theorem 1 above means that typical (as measured by the Standard volume) high-dimensional bipartite 
states are entangled. This statement is of global nature, and is therefore not helpful when one faces a 
given state and has to decide about its separability. In turn, lower bounds for the inradii of sets E, or 
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for the quantity e(7í) mentioned above provide a sufficient condition for separability which, while very 
useful in some contexts, is conclusive for only a very small fraction of separable states. A very simple 
and efficient (at least in small dimensions) necessary condition is the so-called Peres partial transpose 
criterion ([3]). We briefly recali what partial transposition is, say when the state space is H = C D <g> C D . 
Fix a basis (ei, . . . ,eu) of C D , for example the canonical basis. Any state p on Tí can be uniquely 
expressed as 

D 

P= 2J Pia,jp\e-i ® e a )(ej <g> ep\ 

i,j,a,@= 1 

We then consider Tp, the partial transpose of p (with respect to the first subsystem), dcfined by the 
formula 

{Tp) ia ,jP = Pja,i/3 

One checks that Tp is also an Hermitian operator of trace one which, however, is not necessarily positive. 
Following we write PPT (Positive Partial Transpose) for the set of states p for which Tp is also 
positive. Note that while Tp depends on the choice of the basis, its eigenvalues do not; so the set PPT is 
basis-independent. The Peres criterion simply reads: any separable state has a positive partial transpose. 
In other words, we have the inclusions S C PPT C T>. An natural question (asked for example in [2]) 
is then: what is the volume of PPT? Is it close to the volume of S, or rather to the volume of T>1 The 
following thcorcm answers the latter question. 

Theorem 4 (Asymptotic roughness of PPT criterion) There exists an absolute constant cq > 
such that for any bipartite system {Ti = C D <g> C D ) 

An immediate corollary to Theorems 1 and 4 is that, for large D, the volume of PPT states is much 
bigger than the volume of separable states. This means that a typical high-dimensional PPT state is 
entangled, and thus that Peres criterion becomes less and less precise as the dimension increases. Theorem 
4 is just a sample result; similar analysis may be performed for other bipartite and multipartite systems. 
On the other hand, it is quite conceivablc that the quantity from Theorem 4 admits a nontrivial (i.e., 
< 1, independent of D) uppcr bound, even in the case of bipartite systems. [This would imply that vast 
majority of entangled states are detected by the PPT criterion.] 

2 Preliminaries about convex bodies 

In this section we list some known (or elementary) facts from convexity theory that will be needed in the 
sequcl. A convex body is a convex compact subset of R™ (or of C m identified with R 2m , or of an affme 
subspace of R") with nonempty interior. We will say that a convex body K is symmetric if it is symmetric 
with respect to the origin (i.e., K = —K). To a convex body K C R™ containing the origin in its interior 
one can canonically associate \\-\\k, the gauge ofK, by setting, for x £ R", \\x\\k '■= inf{A > : x <G XK}. 
If K is symmetric, || • \\k is the norm for which K is the unit ball. If K and L are two convex bodies in 
R™, each of which contains the origin in its interior, we define the geomètric distance between K and L 
as the product a/3, where a and (3 are the smallest numbers for which (3~ 1 L C K C aL . The concept 
generalizcs to bodies with a common interior point contained in an affine subspace of R n ; it quantifics 
how the shapes of K and L differ. 

2.1 The Urysohn inequality 

Let K C R™ be a convex body and u <G 5 n_1 be a unit vector. The width of K in the direction u is 

w(K,u) :=max(u,i) 
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Geometrically, w(K, u) is (up to the sign) the distance from the origin to the hyperplane tangent to K 
in the direction u. The mean width of K is simply the mean of the widths in all directions: 



w(K) := / w(K,u)do~(u) 
Js"- 1 

where a is the Lebesgue measure on the sphere, normalized so that a(S n ^ 1 ) = 1. The classical Urysohn 
inequality, which is of isoperimetric flavor, states that among all convex bodies of given volume, the 
Euclidean ball has the minimal mean width. It can be written as follows (see, e.g., Appendix A of jSJ for 
a proof) 



vol(A') 



where B% is the n-dimensional Euclidean unit ball. The quantity on the left-hand side equals the radius 
of the Euclidean ball which has the same volume as K and, as we mcntioned earlier, is sometimes called 
the volume radius of K; we will denote it vrad(A). We will use the Urysohn inequality to estimate, among 
others, the volume radius of S, which is the quantity we are interested in. It is often more convenient 
to rcwritc the inequality above by expressing the mean width via a Gaussian integral (rather than an 
integral over the sphere): 

vrad(A') < — Emax(i,G), 

where G is aJV(0, Id n ) Gaussian vector, E is the expectation with respect to G and 7„ := \/2T( IL· ^-) /r(^) 
is a constant which is closc to y/n. The identity w(K) = 7^ 1 E max^gx (x, G) is checked by passing to 
polar coordinates. We will occasionally call Emax ie ]f(i, G) the Gaussian mean width of K and denote 
it u>g(K). 

A special situation is when the body K is a polytope with v vertices on the unit sphere. In this case 
it is wcll-known (see Proposition 1.1.3. in ^2] for a nice proof) that the expectation wg{K) is smaller 
than \2\av. A somewhat more complicated argument based on another proof Lemme 0.6.2) and 

supplemented by numèrics shows that the same estimate works for a symmetric polytope with 2v vertices 
on the unit sphere (or just in the unit ball) provided v > 1. We eventually obtain that, for v > 1 points 
{ x i)i=i m the unit ball, 

vïa,d(conv{±Xi}" =1 ) ^ w(conv{±xJ^ =1 ) < — V2lnv. (2) 

7n 

Note that since j n = wqÍB^) ~ \Jn 1 the quantity on the right will be small as long as lnu <C n. 



2.2 Symmetrization and the Rogers-Shephard inequality 

Let H be an affine hyperplane in R™ +1 , not containing the origin, and let W be a convex body in H. 
Consider the symmetrization of W, defined by íl := conv(VT U — W). This procedure is quite natural 
in the framework of convex bodies, since much more information is available about centrally symmetric 
convex bodies as opposed to general ones. In particular, is the unit ball with respect to the norm ||.||n- 
If h > is the distance between H and the origin, then we have the following inequalities 

2 n 

2hvol(W) < vol(O) < 2h vol(W). (3) 

n + 1 

The left hand inequality is an immediate conscqucncc of the Brunn-Minkowski inequality f |14p. and is 
an equality if the body W is centrally symmetric. The right hand inequality is the Rogers-Shephard 
inequality (^H]) and equality is achieved when W is a simplex. The factor ^-y may appear large, but we 
point out that its ?ith root, which is more relevant here, is smaller than 2. 

We now consider symmetric versions of V and S, namely A := conv(2?U —T>) and S := conv(iSU —S) 
(the latter was defined already in the Introduction); then n = d 2 — 1. The set A is just the Hermitian 
part of the unit ball of the usual Schatten class of trace class operators. In other words A = {A E 
B sa {'H) 1 ^ 1}, where ||A||i is the trace class norm, i.e., the sum of singular vàlues of A or, in our 
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Hermitian setting, of the absolute vàlues of the eigenvalues of A. [See section 12.41 for the functional- 
analytic interpretation of £.] The inequalities imply now 



1 / vol(S) \ 1/n / vol(5) \ 1/n / vol(S) \ 1/n 

2^voI(A)J ^ [voKD)J ^ 2 Vvol(A)J (4) 

Therefore, the loss in precision when passing from the ratio of the volume radii of S and T> to the ratio 
of the volume radii of £ and A is at most a factor of 2. To be precise, the exponent obtained for the 
symmetrized bodies is 1/n whereas it should be l/d 2 (which is the reciprocal of the dimension of A and 
£), but this is not at all an issue. First, we can a posteriori change the exponent, this will at most 
slightly affect the constants. Second, the extra n + 1 in © ensures that no modification of the constants 
is actually needed (see 0, Appendix C, for a precise general statement in this direction). 

For future reference we recali here the estimates from [3] for the volume radius of T> (cited in the 
Introduction), for the volume radius of A, and for the mean widths of these sets. 

^«™ d < a ' = (^f))" d 31 (5) 

1 2 

s$ vrad(X>) w(D) < -= (6) 



2Vd Vd 

While, as we mentioned in the Introduction, both a closed expression for vrad(2?) and its precise asymp- 
totic behavior (as the dimension increases) are known, and while presumably similar information about 
vrad(A) may be obtained via the methods of the compact and transparent inequalities from Q and 
© are sufHcient for our asymptotic results. Moreover, having upper bounds for the larger parameter 
w(-) will provide us with additional flexibility. 



2.3 Ellipsoids associated to convex bodies 

An ellipsoid is an affine image of the unit ball B%- There is a one-to-one correspondence between 
symmetric ellipsoids and sealar produets (or Euclidean structures). Indeed, if £ = TB%, with T a linear 
map, thcn £ is the unit ball associated to the sealar product (x,y)s '■= (T~ 1 x,T~ 1 y). If K is a convex 
body in R™, its polar body, denoted K°, is defined by K° = {x £ M™, Vy £ K, (x, y) < 1}. If belongs 
to the interior of K, K° is also a convex body (if not. it fails to be boundcd). We point out that 
||x||ifo = maxj, 6 iï-(a;, y) and that the restriction of the gauge |j ■ \\k° to the sphcrc coincides with w(K, •), 
the width function of K defined in section 12.11 We further note that the polar of a symmetric ellipsoid 
is again a symmetric ellipsoid. The following definition is uscful for our purposes. 

Definition: If K is a convex body in M n , a John resolution of identity associated to K is a finite 
family (#,, Cí)í 6 j, where x% belong to dK n S 1 ™ -1 (if K C B% or B% C K, these are contact points of K 
with the sphcre) and Ci are positive numbers, such that 

1. J2 c i x i = 

2. Vy£m n ,y = J2c i {x i ,y}x i 

Condition 2 can be rephrased using Dirac notation as ^Cj|a;j)(a;j| = Id„. Taking tracc both sides, we 
see that necessarily °i — n - We recali the following result due to John 

Proposition 1 (John's theorem) Let K be a convex body in W l . Then there exists a unique ellipsoid of 
minimal volume containing K, called Ldwner ellipsoid, which we denote Lòw(if). Similarly, there exists 
a unique ellipsoid of maximal volume contained in K, called John ellipsoid, which we denote John(ií). 
Moreover, the equality Lòw(A') = B% holds if and only if K C B% and there exists a John resolution of 
identity associated to K . Similarly, the equality John(A') = B% holds if and only if B% C K and there 
exists a John resolution of identity associated to K . 



6 



For a proof of this theorem, see ^üj. These ellipsoids are affine invariants: for any affine map T, wc 
have John(Tlv) = TJohn(A') and Lòw(TK) = TLòw(K). In the symmetric case (the one we mostly deal 
with in the sequel), John and Lòwner ellipsoids are also symmetric and dual to each other with respect 
to polarity. More precisely, if K is symmetric, we have io\m{K°) = Lòw(K)° and Lòw(K°) = John(-K")°; 
this is a consequence of the fact that for any symmetric ellipsoid £ , we have vol(£°) = vol(i?2 ) 2 /vol(£). 
It is therefore immediate to pass from the statement about John ellipsoid to the statement about Lòwner 
ellipsoid (actually the theorem is usually stated for the John ellipsoid only). 

2.4 Tensor products of convex bodies 

If K and K' are convex bodies, respectively in R n and R™ , their projective tensor product is the convex 
body KèK' in R" (g> R n ' ~ R™"' defined as follows 

K®K' = conv{a; g) x' , x G K, x' G K'}. 

This tcrminology is motivated by the fact that when K and K' are unit balls with respect to some norms, 
the set K®K' is the unit ball of the corresponding projective tensor product norm on R™ (g) R" . The 
relevance of the notion to our discussion is obvious: the set of separable states is the projective tensor 
product of sets of states on factor spaces; more precisely, if H = Hi <8> H2, then S(H) = r D(Hi)é) r D(TÍ2)- 
It is easy to sec that the operation é> commutes with symmetrizations: if K = conv(A' U — K) and 
K' = conv(7T U -K'), then com(K ®K', -K®K') = K®K' . It follows that, in our notation, = 
A(7íi)(Í3A(7Í2)· Sincc, as was pointed out in section l2~2l A is the unit ball in the trace class norm, S(7Y) 
is the unit ball of the projective tensor product of the spaces B sa (l~tk); k = 1,2, each endowed with the 
trace class norm. 

The definition of K®K' - similarly as other definitions, comments and lemmas of this section - 
immcdiately generalizcs to tensor products of any finite number of factors. Howcver, for the sake of 
transpareney, we shall concentrate in this section on the case of two bodies. We also point out that while 
the definition appears to be well-adapted to O-symmetric sets, cones and linear maps as morphisms (this 
is used in the proof of Lemma ^below and, later, in the first paragraph of section 5), the projective 
tensor product is not invariant under afhne maps. We refer to |17| . Chapter 2, for a discussion of related 
categorical issues and to JB] for exhaustive treatment of tensor products of normed spaces. 

If £ = SB2 and £' = S'B% are two ellipsoids, respectively in R" and R n , their Hilbertian tensor 
product is the ellipsoid £ <g>2 £' := {S ® S')B 2 m ■ This definition docs not depend on the choice of S 
and S". It turns out that Lòwner ellipsoids behave well with respect to projective tensor product, as the 
following lemma shows (note that the analogous statement does not hold for the John ellipsoid). 

Lemma 1 Let K C R n and K' C R™ be two convex bodies. Then the Lòwner ellipsoid of their projective 
tensor product is the Hilbertian tensor product of the respective Lòwner ellipsoids. 
In terms of sealar products, for every x,y in R™ and x' , y' in R™ , we have 

{x® x 1 ,y ® y\ àvi{K§)KI) = {x,y) hòw{K) (x' ,y') hòw(KI) 

Proof: First suppose that Lòw(íT) = B 2 and Lòw(íí / ) = B 2 ■ By John's theorem, there exist John 
resolutions of identity for K and K', respectively (xí, Ci)i^i<^k and (x^-, c'^i^j^k' ■ We easily check that 
K®K' C B 2 n = B 2 ®2 B 2 . Moreover, (xí ® x' 3 : , Cic'^i^i^k.i^j^k' is also a John resolution of identity 
for K®K' . It is cnough to check this on pure tensors: take y G R ra and y' G R" and verify the two 
conditions 

!• Ei=i E*=i Cic^í ® x 'j = (Ei=i ciXi) ® (Ej=i Jj x 'j) = ° 

2 - Et=i Ej=i c l c' :j {x l ®x' ó ,y® y')xi ® x' } = (^=i c t {x t , y)x t ) ® (E*=i ^^'j^y')^]) = V®y' 

For the general case, let T and T' be linear maps such that TLòw(K) = B 2 and T'Lòw(A'') = B 2 . 
Using the elementary identities Low(Tüf) = TLòw(K) and (T ®T')(K®K') = (TK)è(T'K'), the result 
follows from the previous particular case. ■ 
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The next "Chevet-Gordon type" lemma relates the mean widths of convex sets to that of their 
projective tensor product. It is most conveniently stated for the mean Gaussian width %(■) defmed in 
section l^n 

Lemma 2 Let K C B% C R" and K' C i?2 C M n be convex sets, one of which is the convex hull of a 
subset of the corresponding unit sphere. Then 

w G {K®K') ïC w G {K) + w G {K'). 

Restating the assertion of the Lemma in terms of the Standard mean width gives 

w{K®K ) ^ w(K) -i w(K ) ^ 



where the second inequality follows from the fact that the sequence (-jk/Vk) is increasing (which can 
be deduced from incqualities proved in |19| : the fact that the cocmcients on the right hand side are 
approximately \j\fn! and l/y/n follows from the simpler relationship "fk ~ Vk.) We thus get seemingly 
strong upper bounds for mean widths and, a fortiori, for volume radii of projective tensor products 
of convex bodies. However, in spite of its elegance, Lemma [21 will play only rather limited role in our 
arguments. This is because, when iterated, it yields asymptotically worse dependence on the number 
of factors than the techniques employed in subsequent sections. Accordingly, its main use will be in 
considerations involving few factors and in ameliorating the numcrical constants which appear in the 
statements of the Theorems. The Lemma can be shown by an argument very similar to that of section 
2.3 of [201, w e skip the proof; here we just note that we do not know whether the hypothesis that one of 
the sets be spanned by a subset of the sphere is necessary. This rcquiremcnt is satisficd for all sets in our 
applications, and it is not at all rcstrictive if the widths of K, K' are estimated via cardinalities of their 
spanning sets. 

Finally, let us note that if K and K' are circled convex bodies in respective complex vector spaces (a 
convex body K C C™ is said to be circled if K = éf^K for all real ■&), then K®K' is circled; as was the 
casc for symmetrizations, the operation absconv (the absolute convex hull) commutes with the projective 
tensor product. 



3 Proof of Theorem 1: Small number of large subsystems 

The upper bound can be obtained through a Standard discrctization argument. Recali that, for 6 > 0, 
a (5-net of a set K is a subset Àí C K such that for each x € K there exists y G N whosc distancc to 
x is smallcr than S. If TV C C D , we write P(AÍ) for the polytope conv{±\x)(x\} x< z^f C B sa (C D ). The 
following elementary lemma shows that if J\í is a net of the unit sphere of C D - which may be identificd 
with S 2D ~ l as a mètric space - then P(AÍ) is a good approximation of A(C D ). 

Lemma 3 Let 6 < fa 0.765 and let TV be a 5-net of the unit sphere of C D . Then 

(1 - 2S 2 + S 4 /2)A(C D ) c P(ÀÍ) C A(C D ) 

Proof: The second inclusion is trivial. Let us check the first one through the corresponding dual 
(polar) norms 

IHIa(CD)° = max (A,±\x){x\)hs = max | (x\A\x) \ = max | A| = || A|| op , 

WMpw = ma ïï\(y\ A \y)\ 

where || • || op is the operator norm and tr(·) is the spectrum (recali that we consider only Hermitian matrices 
here). We need to show that ||A|| PC ^)o > (í-2S 2 +S i /2)\\A\\ op . To this end, let A e B sa (C D ) be such that 
\\A\\ op and the largest eigenvalue of A are both equal to 1, and let x £ C D be a norm one vector such that 
Ax = x. Choose j/o S N verifying ||.t — yo\\ ^ S. We claim that (yo|^4|yo) > 1 — 2<5 2 + <5 4 /2; the inequality 
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between the norms follows then by homogeneity. The claim is easily established by writing y — y\ + y 2 
with Vl = (y Q \x)x and noting that ( Vl \A\yi) > \(y \x)\ 2 > (l-S 2 /2) 2 , \(y 2 \A\y 2 }\ < \\y 2 \\ 2 < 'ó 2 (l - <5 2 /4) 
tmd{ yi \A\y 2 )=0. M 

Tensoring the conclusion of the preceding lemma and recalling that A(C D )® N = S(7í) yields an 
inclusion (1 -2S 2 +S 4 /2) N E(H) C P{M) èN . Observe that P(AÍ) èN is anothcr symmctric polytopc with 
at most 2(#AT) N vertices. It is wcll-known ([21], Lemma 4.10) that for S ^ 1 we can find ó-nets in S 2D ~ l 
of cardinality not exceeding ^ (1 + 2/8) 2D . We thus get a bound on the volume of the polytope using 
the estimation 

We note for future reference that the above discussion can be carried out using other Euclidean 
structures, the only constraint being that the vertices of the polytope are of norm not exceeding 1. Thus, 
if £ is any ellipsoid containing £(7i), we obtain 



ol(S(7í)) \ 1/rf2 
vol(f) ) 



inf 



0< 5 <V^7^(l-2<5 2 + íV2) 



^21n((l + 2/í) 



2DN\ 



(7) 



Using £ = B HS and, say 6 = l/ViVln2iV lcads to an estimate vrad(E(H)) = 0((DN In N) l l 2 /d). 
[Indeed, for large d we have then 7d 2 - d and, for large TV, (1 - 2S 2 + 8 A /2) N ~ 1 and ln ((1 + 2/5) 2DN ) ~ 
DNhiN.] We obtain the upper bound announced in Theorem 1 by combining this with the lower bound 
on vrad(A(7í)) given by JSJ) and with Finally, the more precise statements about the constants 
indicatcd in the introduction can be obtained by using JJJ, Lemma [21 and Appendix C of PJ- 

The lower estimate will be proved by showing that S(7í) contains a Euclidean ball of appropriately 
large volume. We start with the following lemma 

Lemma 4 Let K := (B®)®" 1 be the projective tensor product ofm copies of the D-dimensional Euclidean 
ball (real or complex). Then K contains the following múltiple of the Euclidean ball in (C D )® m (resp., 
(JH D )® m ) 

D (m-1)/ 2 B ? CK - 

Proof: Let (ei,...,en) be the canonical basis of C D . We write an element A of (C- )®" 1 as a 
generalized matrix: A = (cii 1 ...i m ) stands for 



A 



E 



'1 ! 



= 1 



The norm dual to 



K IS 



\A\\ko 



max 
x\...,x m €B? 



E 



Let us write ||A||2 for the Euclidean norm of A, i.e., \\A\\\ = ^2 



(8) 



I 2 . We want to show that 



VA G (C 15 )®" 1 we have D _ ( m_1 )/ 2 m|| 2 s$ \\A\\ K o. By homogeneity, we can assume that \\A\\ 2 = 1. Now 



choose x 



randomly and indcpcndently to be uniformly distributed on the unit spherc of 



We write E for the corresponding expectation. Let X^, k = 1, . . . , D, be the random variable obtained 
by summing only on the fcth "hypcr-slice," that is 



X k =X k {x\ 



l ):= 



D 

E 

, im-l = 1 



jm— 1 



We can easily calculate E|Xfc| 2 since many cancellations come from the fact that E(xfx' ) = j)S k i6ij. We 
obtain 

1 D 

/ , | a ii...i„_ifc| 



nx k \ 2 



D' 



,i m -i = 1 
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Therefore, E(^^ =1 |Xj,| 2 ) = . This implies the existence of unit vectors x 1 , . . . ,x m 1 such that, 

denoting Yfc := Xk(x l , . . . , x m ~ l ), we have X)fc=i l^fe| 2 ^ d™- 1 ■ Choosing these points in JÏÏJ), we obtain 



1 /2 

2 k=l \k=l / 



Remark: It turns out that this lcmma is surprisingly sharp. As shown in Lemma IHlbclow, even the a 

priori larger volume radius of K is (up to multiplicativc factor depending only on m) of the same order 
D -( m -i)/2_ Note that for 

m = 2 we obtain even the optimal constant (these are the classical inclusion 
relations for balls in Schatten classes). 

Let Tí = (C £> )' 8lJV and T(TL) C B(TL) be the convex hull of rank one product operators 
Y{H) ■.= com{\xi®···®x N ){y- í O···Oj/jvI : xi,yi,-- • ,xn,Vn G B$ }. 

The convex body T(TL) is most naturally seen as the Nih projective tensor power of the set of (not 
necessarily Hermitian) operators on C D whose trace class norm is < 1. However, it can also be identificd 
with (i?! 3 )® 2 ^ when we identify B(fC) with K. ® K.. It then follows from the preceding lemma that T(TL) 
contains a Euclidean ball (a Hilbert-Schmidt ball in this context) of radius £) _ ( 2JV_1 )/ 2 . 

The next lemma will relate T(TÍ) to S(7í), allowing to dcduce that the latter body contains a suitably 
large Euclidean ball. Let be the geomètric distance between the sets A and £ corresponding to N 
qubits, i.c. the smallest positive number such that A((C 2 )® Ar ) C d N Y,((C 2 )® N ). We then have 

Lemma 5 Let n : B(TÍ) — > B sa (TL) be the projection onto Hermitian part, ir(A) := \{A + A'). Then 

n(T(H)) C d N Z(H). 

Proof: It is enough to show that extreme points of 7r(r) are contained in d^T,. Any extreme point 
A of 7r(r) can be written as 

A = - (\x! ® • • • ® x N )(yi ® • • • ® y N \ + \y x ® ■ ■ ■ ®y N )(xi ® • • • ® x N \) 

It may appear at the first sight that the above representation shows that A is separable. However, whilc 
the two terms in the parentheses are indeed product operators, they are not Hermitian and we can only 
conclude that A e A(W) (as a Hermitian operator whose trace class norm is < 1). 

Let Tíi be the 2-dimcnsional subspace of C D spanncd by Xi and yi (if the vectors are proportional, 
add any vector to get a 2-dimensional space) and let H' := Hi ® ■ • • ® Wjv- Then A can be considered 
as an operator on Tí'; more precisely, as an element of A(TÍ') (and, conversely, any operator acting on 
H' can be canonically liftcd to one acting on Tí). Note that since we are in the asymptotics wherc N is 
small and D is large, the dimension of the problem has been dramatically decreased. Sincc A belongs to 
A(H'), it also belongs to dN^{TÍ'), and thus to c?atE(7í). ■ 

To effectively apply Lemma|Sl upper bounds on d^ are needed. It has been observed in [3], Appendix 
H, that dN < 6 W / 2 . This estimate has been subsequently slightly improved in 0] to 2/3 x 6 W / 2 , and 
even to ajy x 6 N ^ 2 with liniAr-foo oat w 0.624. [These bounds are in fact based on estimates for the 
a priori larger quantities, the geomètric distances between E((C 2 ) lglAr ) (or S((C 2 )® N ) — ld/2 N ) and 
Hilbert-Schmidt balls, which are essentially reciprocals of the Hilbert-Schmidt inradii of these sets; and 
in this latter context they are "nearly optimal;" see Theorem 3 and the comments following it.] Thus 
S(7í) D 3/2 x 6~ Ar ^ 2 7r(r(7í)); combining this inclusion with Lemma 01 or, more precisely, with the 
observation following the definition of T(7í), we conclude that E(7í) contains the Hilbert-Schmidt ball of 
radius 3/2 x q-n/2 D -(2N-i)/2 = 3 / 2 x Q-N/2/ d i-i/2N_ This givcs a lower bound for the volume radius 

which, together with JSJ, yields 

vol(E) V /d2 > 3/4 x 6^/ 2 



vol(A) ) " dV2-l/2JV 
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To conclude the proof of the theorem we pass to the bodies S and T> using 

Remarks: (1) The bipartitc casc (N = 2) was studied in JT], which contains in particular the 
remarkable result that D and S have then the same inradius, which equals 1/ y^d(d — 1). However, the 
inradii of the symmctrizcd bodies, which arc (in both cases) comparable with the volumc radii, are 
(asymptotically) very diffcrent: the inradius of A equals 1/Vd, while the inradius of E has just been 
shown to be of order l/d 3 ' 4 . 

(2) The upper bound on cIn was further slightly improved in subsequent versions of H] and, most 

recently, in [22], where it was shown that S contains a ball of radius \fï x Q~ N / 2 , which implies Ún < 

( 3 JV-i (2 ív _ 1)) i/2_ 

4 Proof of Theorem 2: Large number of small subsystems 

This part deals with asymptotic estimations of volumes when N, the number of subsystems, tends to 
infinity whereas the dimension D of each subsystem remains bounded. This case is probably the one with 
most physical interest. When all subsystems are 2-dimensional (qubits), Theorem 2 has been proved in 

El 

For the proof of the present case it is more convenient to deal with an affine image of £ which is 
more balanced than £ itself or, equivalently, to consider a new sealar product which is better adapted to 
the analysis of £ than the Hilbert-Schmidt product. It is a known phenomenon that high-dimensional 
convex bodies often enjoy hidden symmetries (or approximatc symmetries) that are only revealed when 
one looks at them through a suitablc Euclidcan structurc. In our situation, in the asymptotics when 
N is large and D is bounded, our Theorem 1, obtained with the usual Euclidean structure, misses the 
genuine volumic behavior of the bodies S and £. The appropriate sealar product will be here the one 
derived from the Lòwner ellipsoid of £. To see why, just look back at the equation Q. This bound is 
the tightest when £ is the Lòwner ellipsoid of £. It turns out that it is possible to describe completely 
this ellipsoid. First recali that £(7í) = A(C D )&> . . . <è)A(C D ). By Lemma^] it is enough to calculate the 
Lòwner ellipsoid of A(C D ), and by tensoring we get the Lòwner ellipsoid of £(W). 

We will determine the Lòwner ellipsoid of A(C £> ) using the following elementary general lemma. 

Lemma 6 Let u be a unit vector in R n and Hq = {u}- 1 be the orthogonal hyperplane. Let h > be real 
number and H = Hq + hu be a translate of Hq. We consider a convex body fi = cotïv{W U — W}, where 
W is a convex body in H . Then the following assertions are equivalent: 

1. íl is in Lòwner position (i.e. the unit ball B% is the Lòwner ellipsoid ofíï). 

2. h = l/\/n and i?2 l~l H is the Lòwner ellipsoid ofW. 

Proof: Assume that Í2 is in Lòwner position, so fl C B%. By John's theorem (Proposition^l, this 
means that there exists a John resolution of identity (xí,Cí), with xí G n S 71 " 1 and ^Cí\Xí)(xí\ = Id. 
Since Xí are extreme points of íi, thcy must be in W or — W , and since is symmetric we can assume they 
are all in W . Let P be the orthogonal projection onto Hq. It follows that ^ Cí\Pxí)(Pxí\ = P; note that 
ll-Pïill = Vi — h 2 . This mcans that (Pxí, c') is a John resolution of identity in Hq, with cj ; = (1 — h 2 )ci. 
By PropositionQ (1 - h 2 )~^ 2 PW is in Lòwner position, or equivalently H H is the Lòwner ellipsoid 
of W. Moreover, we have the conditions ^2 d i — n — 1 and "Y^a = n, which force h = 1/^/n. 

The converse follows by retracing the above argument in the oppositc direction. ■ 

Lemma 7 The Lòwner ellipsoid of A(C D ) is determined by the following sealar product defined for 
A,BeB sa (C D ) by 

<A£0 L5w(A(CD)) = (1 + -p)tr(AB) - ^tr(A)tr(B) 

Proof: First remark that, because of its uniqueness, the Lòwner ellipsoid of any body K must 
inherit all the symmetries of K. Since Ï?(C £> ) is invariant under unitary conjugations (i.e., the maps 
ivu{X) = UXW with U G U(D)), the ellipsoid Lòw(D'C D )) must be similarly invariant, and thus it is 
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the circumscribed Euclidean ball. This is actually somewhat delicate; the relevant fact is that the action 
of the group {ttjj : U £ U(D)} on the hyperplane T\ containing V(C D ) is irreducible (in particular, Id/D 
is the only fixed point and so it must be the center). In the language of convex bodies (see [2H]) we say 
that T>(C D ) has enough symmetries. On the other hand, the fact that this group acts transitively on the 
set of extreme points of T)(C D ) (which are the pure states) implies by itself only that all these extreme 
points are contact points, i.e., belong to the boundary of the Lòwner ellipsoid. 

By Lemma El we know that the trace of the Lòwner ellipsoid of A(C D ) on the hyperplane T\ is a 
Euclidean ball (centered at Id/D), so we have for some real numbers a and (3 

(A, B) Low(A{C o }) = atr(AB) + í3tr(A)tr(B) 

To determine a and (3, observe that (/?, /?)low(A(c d )) = 1 f° r am/ pure state (hence contact point) p, and 
that (from the condition on h in Lcmnia[üJ note that n = D 2 ), (Id/D, H/D)low(A(c- d )) = 1/D 2 . This 
gives a=l + l/D and /3=-l/D. U 

It is possible to use the Lòwner ellipsoid of as the ellipsoid £ in formula J7J). Choosing 6 = 

Í/VN \n2N as earlier, we obtain the upper estimatc in the formula 

1 ^ / vol(S) \ 1/d2 ^ CVDNlnN 



d Vvol(Lòw(S)) 

The lower estimatc follows from the classical fact (JH) that for a centrally symmetric convex body K 
in R™, we have the inclusion Lòw(K) C y/n K. Improvements on constants can be derived from a result of 
Barthe ([23) which asserts that for an n-dimensional symmetric body K, the ratio (vol(K )/vo\(Lòw(K)) 1 / n 
is minimal when K is the unit ball of P{ (and thus is at least (4/ 7 r) 1 / 2 (r(n/2 + l)/T(n + l)) 1/n ) = 

^J^{l - i we reca11 that here n = d2 )- 

It remains to calculate the volume of the Lòwner ellipsoid of S. It follows from Lcmmanjthat, if $ 
denotes a linear map which sends B% onto Lòw(A(C £> )), then := maps Bhs onto Low(E(W)). 
By Lemma we can define $ by ®(A) = (1 + l/D)- l l 2 A if A has trace and $(Id) = \f~Dlà. This 
gives det($) = ^ D{1 + l/D) 1 - ' 2 and finally 

VOl( ^ Ò D W(E)) = det(*) = detW^ 2 "- 2 = d-^ 2 (10) 
vol(B HS ) 

To complete the proof it remains to combine the volume estimates (|5|).<|üll. (|10|) and 

5 Proof of Theorem 3: Tight upper bounds on the inradii 

The idea is to show that there is a projection of the set S = S((C 2 )' 8JV ) which is dcmonstrably small. 
Spccifically, lct P : B sa (C 2 ) — > B sa (C 2 ) be the orthogonal projection onto the 3-dimensional subspace of 
trace zero matrices and let l·l := P® N . Since projective tensor products commute with linear maps (in 
the sense that (Ti ® . . . ® T a )(iï"i(8) . . . ®K„) = (Tiifi)® . . . <g)(T s K s )), we have 

n (^((C 2 )^)) = P® N (A(C 2 )^) = (P(A(C 2 ))f N . 

Sincc A(C 2 ) is a cylindcr whose axis coincides with the kerncl of P, the imagc of A(C 2 ) under P is a 
3-dimensional Euclidean ball of radius 1/V2 (more precisely, V(C 2 ) — Id/2). Accordingly, the image of the 
corresponding set £ under n is the nth projective tensor power of such a ball. The relevant parameters of 
such sets (mean widths, volume radii) are majorized by the following statement, which is a complement 
of Lemma 4. 

Lemma 8 There is an absolute constant C\ such that if D, m are positive integers (m ^ 2) and K := 
(_B,P)® m is the projective tensor product ofm copies of the D-dimensional Euclidean ball (real or complex), 
then 

,,-rrs , t ^n C\\Jm ln m 
w3d(K)^w(K)< ^ m _ 1)/2 
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Proof: The first inequality in Lemma[S]- comparing the volume radius of a convex body to its mean 
width - is just the Urysohn inequality. The proof of the second follows the same lines as - but is simpler 
than - the approach that yielded formula (JJJ . We present the details only for the case that is relevant to 
Theorem 3, namely D = 3 in the real setting. The argument is very similar to the one used in [3] to obtain 
an upper bound for the volume of E((C 2 )® m ), where the reader is referred for additional background and 
justifications. 

We start by noting that, for any 5 £ (0, y/2), the sphere S 2 contains a (5-nct T whosc cardinality is 
^ 16/(5 . On the other hand, an element ary geomètric argument shows that the convex hull of any <5-net 
contains a (Euclidean) ball of radius 1 — <5 2 /2. Accordingly, the convex hull of J : ® m , the mth tensor 
powcr of T, contains (1 - <5 2 /2) m (Bf )® m . Sincc #(T® m ) s$ (#.F) m < (4/<5) 2m , it follows from © that 

w{{1 - 5 2 /2) m (Bïf m ) ^ w{conv{F® m )) < 73^\/21n((4/(5) 2m ) 

and, consequently, 

vrad((Bl) ém ) < w{(B%)® m ) < 2^(1 - 5 2 /2)- m ^/m\n(4/S). 

We now choose 6 = \j\frn ln2m; then, for large m, (1 — 5 2 /2)~ m ~ 1, ln (4/(5) ~ 4 lnm and 73™ ~ 3 m / 2 . 
This yields the asserted upper bound on the mean width with C\ ~ \/2/3 for large m. A more careful 
calculation using Lemma El for small m shows that C\ = 1.673 works for all m. [The discussion of 
constants above and in the remainder of this section is specific to D = 2, but can be analogously carried 
over in the general case; for examplc, it is easy to show that C\ ~ 1 works for large m and all D] As 
in PJ, further improvements arc possible by working with explicit nets and/or by noting that T can be 
chosen to be symmctric, but we will not pursue them since it is likcly that, for the problcm at hand, 
sharper estimates can be obtained by analytic mcthods. ■ 

We now return to the proof of Theorem 3. From prior considerations, the image of E under the 
projection II is congruent to 2~ N l 2 (B\)® N and so, by Lemma|Hl 

w(ILÇE)) = 2- N / 2 w((Bïf N < 3 1 / 2 C 1 VMnlV 6-^. 

Since the inradius is trivially smaller than the mean width, and since the inradius of the set S of separable 
states does not exceed that of its symmetrization E, Theorem 3 follows (with Cq = 3 X I 2 C\ < 3). It is 
also routine to obtain an (identical or ncarly identical) upper bound on the mean width of S; however, 
some care is needed because the dimensions of S and E are not the same. 

Since our argument is of a global nature, we are not able to produce an explicit state in the ball 
of radius r := Co V N ln N 6~ N ^ 2 and centered at Id/d which is not separable. However, Standard 
tools of asymptotic geomètric analysis (cf. [7j or [21\ ) imply that even stronger properties are satisfied 
"generically." Denote by TÍq the do := 3^-dimcnsional range of the projection n. Then, for most of the 
unit vectors u on the unit sphere of Tío, the state er = Id/2^ + 2ru is not separable, and the same holds 
for any state p such that II(p) = H(a) (here "most of means that the exceptional set is of normalizcd 
measurc < 2~ d °/ 2 ). Moreover, a slightly weaker but comparable property holds simultaneously for all 
u in a "genèric" (in the sense of the invariant measure on the corresponding Grassmannian) [~g£o/2~|- 
dimensional subspace of 7io ■ 

As mentioned in the introduction, the above argument can be generalized to other vàlues of D and 
leads to an upper bound for the inradii of the sets E and S corresponding to (C D )® N which is, up to a 
logarithmic factor, (D(D + \))~ N / 2 . While this upper bound differs, for D > 2, from the lower bound for 
the inradius of S that was found in 0, it gives the correct order for the inradius of the set E. Indeed, E 
contains Lòw(E)/d, which itself contains {D {D + 1))~ n / 2 Bhs (this argument parallels the one for qubits 
given in Appendix H, where we refer the reader for details). 

6 Proof of Theorem 4: Asymptotic roughness of PPT criterion 

We will use the following proposition which is taken from [25) . 
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Proposition 2 Let K and L be convex bodies in W 1 with the same centroid. Then 

vol(A")vol(i) < vol(K - L)vol{K n L), 
where K — L := {x — y : x € K,y G L} is the Minkowski difference. 

We apply the proposition with K = V and L = TV. This gives 

vpi(x>) vpi(x>-rx>) 

vol(PPT) vol(P) 1 ' 

Note that we actually used the proposition in the space 7| ~ W 1 , seen as a vector space with Id/d as 
origin; the point Id/d is also the centroid of V (and TV), it is even the the only fixed point under the 
group of symmetries of V (cf. the proof of Lcmma[7|>. We rewrite (|11J) using volumc radii 

/ vol(X>) \ V " vrad(r>-TX>) 
Vvol(PPT) ) vrad(P) 

Recali that we majorized the volume radius by the mean width using the Urysohn inequality (here 
mean width is also taken in T\). This is vcry convenient since, as one can check immediately from 
the definition, the mean width is additive with respect to the Minkowski operations: w(V — TV) = 
w(V) + w(—TV). Since — T is an isometry, it preserves the uniform measure on the sphere S*™" 1 and 
thus w(-TV) = w(V). This gives 

/ vol(P) \ V " 2w{V) 
Vvol(PPT) ) ^ vrad(P) 

We now appeal to (JüJ) to majorize the right hand side by 8, which proves Theorem 4 with co = 1/8. 
Note also that one can use the explicit formula Q for the volume of V; this shows that asymptotically 
(the lower estimate for) cq tends to e -1 / 4 /4 w 0.195. Moreover, numerical evidcncc (which presumably 
can be confirmed by rigorous calculation using appropriatc Stirling-typc formulae, or incqualities similar 
to those of ^ü]) suggests that vrad(X>) ^ e" 1 / 4 j\fd for all dimensions d, which would imply that the 
asymptotic value is in fact a bound. 

Finally, observo that the above proof works actually for any (Hilbcrt-Schmidt) isometry in Ti which 
fixes Id/d, and not only for the partial transpose T. 
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